Existence theorems in the geometrically non-linear 
6-parametric theory of elastic plates 

Mircea Birsan* and Patrizio Neff ' 

O 

(N 

!— i ■ March 12, 2013 

a 



oo 



PL, 
< 



X 



Abstract 



In this paper we show the existence of global minimizers for the geometrically 
exact, non-linear equations of elastic plates, in the framework of the general 6- 
■ parametric shell theory. A characteristic feature of this model for shells is the 

appearance of two independent kinematic fields: the translation vector field and the 
rotation tensor field (representing in total 6 independent scalar kinematic variables). 
For isotropic plates, we prove the existence theorem by applying the direct methods 
CN , of the calculus of variations. Then, we generalize our existence result to the case 

of anisotropic plates. We also present a detailed comparison with a previously 
Q\ ■ established Cosserat plate model. 
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1 Introduction 

In this paper we present an existence theorem for the geometrically non-linear equations 
of elastic plates, in the framework of the 6-parametric shell theory. 

The general (6-parametric) non-linear theory of shells, originally proposed by Reissner 
has been considerably developed in the last 30 years. This theory and the most 
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results in the field have been presented in the books of Libai and Simmonds [24] and 
recently Chroscielewski, Makowski and Pietraszkiewicz [TU]. The model is based on a 
dimension-reduction procedure of the three-dimensional formulation of the problem to the 
two-dimensional one, and is expressed through stress resultants and work-averaged defor- 
mation fields defined on the shell base surface. Thus, the local equilibrium equations for 
shells are derived by an exact through-the-thickness integration of the three-dimensional 
independent balance laws for linear momentum and angular momentum. The deforma- 
tion of the shell is characterized by two independent kinematic fields: the translation 
(displacement) vector and the rotation tensor. The appearance of the rotation tensor as 
an independent kinematic field variable is one of the most characteristic features of this 
general shell theory. In this respect, we mention that the kinematical structure of the 
non-linear 6-parametric shell theory is identical to that of the Cosserat shell model (i.e., 
the material surface with a triad of rigidly rotating directors attached to any point) pro- 
posed initially by the Cosserat brothers [H] and developed subsequently by Zhilin [4"2] . 
Zubov [13], Altenbach and Zhilin [3], Eremeyev and Zubov [20], Birsan and Altenbach [7J, 
among others. A related Cosserat shell-model has been establish recently by Neff [271 [29] 
using the so-called derivation approach. 

The subject of derivation and justification of the non-linear, geometrically exact equations 
for plates and shells has been treated in many works, and the existence of solutions has 
been investigated using a variety of methods, such as the method of formal asymptotic 
expansions or the T-convergence analysis, see e.g. [ffil EH [391 E El El ESI [31]. For an 
extensive treatment of this topic, as well as many bibliographic references, we refer to the 
books of Ciarlet [T21 [33] • To the best of our knowledge, one cannot find in the literature 
any existence theorem for the non-linear 6-parametric theory of plates or shells developed 
in [211 110]. In our work, we describe the non-linear equations of elastic plates as a two- 
fields minimization problem of the total potential energy and we prove the existence of 
minimizers by applying the direct methods of the calculus of variations. The first result 
concerning the existence of minimizers for a geometrically exact (6-parametric) Cosserat 
plate model has been presented by the second author in [271 EH]. Due to differences 
in notation, this result has not been much noticed. In this paper, we employ similar 
techniques as in [27] and adapt the existence proof to the general 6-parametric plate 
equations. 

In the framework of the linearized 6-parametric theory, the existence of weak solutions for 
micropolar elastic shells has been proved recently in [TB] . We mention that the kinematic 
structure of the general 6-parametric shell theory differs from that of the so-called Cosserat 
surfaces, i.e. material surfaces with one or more deformable directors [221 EHl H]. In 
particular, the kinematics of Cosserat surfaces with one deformable director [261 [38] leaves 
indefinite the drilling rotation about the director, while the general 6-parametric shell 
theory is able to take into account such drilling rotation. For the linear theory of Cosserat 
surfaces, the existence theorems have been established in [T5l 151 16]. 

Here is the outline of our paper: In Section [2] we briefly review the field equations of 
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the non-linear 6-parametric plate theory. Then, in Section [3] we prove the existence 
theorem for isotropic plates. The generalization of the existence result for anisotropic 
plates is presented in Section HI We also show that the existence theorem remains valid 
in the case of some alternative relaxed boundary conditions for the rotation field. In 
Section [5] we present a detailed comparison between the non-linear 6-parametric plate 
and the Cosserat plate model proposed and investigated by the second author in [27J 
[25] . Although this Cosserat model for plates has been obtained independently by a 
formal dimensional reduction of a finite-strain three-dimensional micropolar model (see 
also [3 1 J ) , we show here that the kinematical variables and the strain measures of the 
two models essentially coincide. Moreover, the expressions of the elastic strain energies 
become identical, provided we make a suitable identification of the constitutive coefficients 
for isotropic plates in the two approaches. 

The linearized version of this model has also been investigated in [23 |2§J [2S1 123 E2] and its 
relations to the Reissner-Mindlin, Kirchhoff-Love and other classical plate models have 
been discussed. 



2 Basic equations of geometrically exact elastic plates 

The governing equations of the general 6-parametric non-linear theory of shells have been 
derived in [2U QUI E] by direct integration of the two independent fundamental principles 
of continuum mechanics: the three-dimensional balance laws of linear momentum and 
angular momentum. In this section we summarize these equations, specialized here for 
the case of plates. To this aim, we employ mainly the notations introduced in [2H (TTJ [10] . 
Consider an elastic plate, which is a three-dimensional body identified in the reference 
(undeformed) configuration with a region D, — {(xi, x 2 , z) \ (xi, x 2 ) G u), z G [ — | , | ] } of 
the Euclidean space. Here h > is the thickness of the plate and u C M 2 is a bounded, 
open domain with Lipschitz boundary duj. Relative to an inertial frame (O, e»), with e» 
orthonormal vectors (i = 1,2,3), the position vector p of any point of Q can be written 
as 

p(x a , z) = x + ze 3 , x = x a e a , (x 1 ,x 2 )euj, ze [--,-]. ^ 

Throughout this paper, we employ the usual convention of summation over repeated 
indices. The Latin indices i,j, ... take the values {1,2,3} and the Greek indices a,/3, ... 
range over the set {1,2}. 

In the deformed configuration, the base surface of the plate (shell) is represented by 
the position vector y = x(x), where x '■ u G M 2 — > IR 3 is the surface deformation 
mapping. Let the vector field u(x) represent the translations (displacements) and the 
proper orthogonal tensor field Q(x) designate the rotations of the shell cross-sections. 
Then the deformed configuration of the plate is given by 

V = X(x) = x + u(x), d % = Qe u i = 1,2,3. (2) 
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Figure 1: The reference base surface u of the plate and the deformed surface 5, described 
by the surface deformation mapping y = x( x ) an d the independent rotation tensor field 
Q(x). 



The vectors di introduced in (T5]) are three ortho normal directors (see Figure 1) attached to 
any point of the deformed base surface S = x( u )- Thus, the rotation field Q(x) G 50(3) 
can be written as 

Q = di®ei. (3) 

According to the Lagrangian description, let / and c be the external surface resultant 
force and couple vectors applied at any point y G 5, but measured per unit area of u. 
Also, let n v = Nu and ra v = Mv be the internal contact stress and couple resultant 
vectors defined at an arbitrary boundary curve dG C 5, but measured per unit length 
of the undeformed boundary dj C ui. (We have denoted here by G = x(l) an d v 
is the external unit normal vector to lying in the plane of u.) Here, the tensors 
N = N ia ei eg) e a and M = M ia ei <g) e a are the internal surface stress resultant and stress 
couple resultant tensors (of the first Piola-Kirchhoff stress tensor type). Then, the local 
equilibrium equations are given in the form [2o"| [TP] 

Div s N + f = 0, Div s M + ax\(NF T - FN T ) + c = 0, (4) 

where F = Grad s y = y a ®e a is the surface gradient of deformation and Div s iV = N ia ^ a ei 
, Div s M = Mi a ^ a ei . As usual, we denote the partial derivative with respect to x a by 
(•),« — af - (')• I n © the superscript ( -) T denotes the transpose and axl(A) is the axial 
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vector of any three-dimensional skew-symmetric tensor A, given by 

axl(A) = A 32 e 1 + A 13 e 2 + A 21 e 3 , for A = <g> e 3 - , A T = -A, (5) 

such that Av = axl(A) x v, for any vector v. The corresponding weak form of the local 
balance equations has been presented in [21] Chap. VIII, or in [T7] . 
To formulate the boundary conditions, we take a disjoint partition of the boundary curve 
du = dojd U dujf , dojd fl duf = 0, with length (dud) > 0. We consider the following 
boundary conditions [TUl ES] 

u — u* = 0, Q — Q* = along <9a>d , (6) 

iW - n* = 0, Mv - m* = along du;/ , (7) 

where n* and m* are the external boundary resultant force and couple vectors applied 
along the part duf of the boundary du. In general, from a modelling point of view, it is 
a difficult task to specify boundary conditions for the rotation Q in (JB} at the boundary 
dua ■ In [T?| the following requirement is considered: in relations ()6]) the functions u* and 
Q* defined on dud should be found from the Dirichlet boundary conditions u 3D (x a , z) = 
ul D (x a , z) for the three-dimensional body fl, at any point (x a , z) e dQd = dud x [— |, |]. 
Thus, the functions it* and Q* should be determined from the condition that the work 
done along dud by the resultant stress and couple vectors n„ , ra v on the translation u* 
and rotation Q* be the same as the work done along dQd by the nominal three-dimensional 
stress vector t u (x a , z) on the translation u\ D {x a , z). 

In the general resultant theory of shells, the strain measures are the strain tensor E and 
the bending tensor K, given by [TUl El EE] 

E = Q T (e a ® e a ), £ a = y a -d a , (8) 

K = Q T {x a <g> e a ), x a = &xl(Q a Q T ). (9) 

We mention that the kinematical structure OH]),© of the general shell theory is identical 
with that of the classical version of the Cosserat shell [El H21 [31 EHl E]. In the case of 
plates, these strain tensors can be written in component form relative to the tensor basis 
{ej Cg) e a } as 

E = E ia ei ®e a = (y a ■ d { - 6 ia )ei <g> e a , (10) 

= (da,« ■ d 3)ei g) e Q + (d 3ja • d 1 )e 2 ®e a + (d ljCe ■ d 2 )e 3 <g>e a , 

where 5i a is the Kronecker symbol and is the permutation symbol. 

According to the hyperelasticity assumption, the constitutive equations for elastic plates 

are given in the form [2H [101 [18] 

dW dW 
N = Q—, M = Q— , (12) 
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where 

W = W(E,K) (13) 

is the strain energy density. Then, the equations (jlj) are the Euler-Lagrange equations 
corresponding to the minimization problem of the total potential energy. In order to 
characterize the material of the elastic plate, one has to specify the expression of the 
potential energy function (1131) . In the paper [TS] the conditions for invariance of the strain 
energy under change of the reference placement are discussed and the local symmetry 
group is established. The structure of the local symmetry group puts some constraints on 
the form of W, which allows one to simplify the expression of W. From this representation, 
the strain energy density corresponding to physically linear isotropic plates is given by 
(see [IS], Sect. 10) 



W(E, K) = W mh {E) + W hcnd (K), 
2W mh (E) = aitr 2 £|| + a 2 trE\ + a 3 tr(EqE\\) + a 4 nEE T n, 
2W hend {K) = /?itr 2 K|, + /3 2 trKj + /3 3 tr(KjK\\) + f3 in KK T n, 
E\\ = E - {n®n)E, K\\ = K - (n <g> n)K , n = e 3 , 

where the coefficients a& , ftk (k — 1,2,3,4) are constant material parameters. 



(14) 



Remark 1. In the works [774 HO]/ the authors have employed a particular form of the 
expression (IT4"|) for the strain energy density in the isotropic homogeneous case, namely 



2W(E,K) = C[vti 2 E\\ + (1 - v) tr{E\ E { \)] + a s C(l - v)nEE T n 

+ D[vti 2 K\\ + (1 - u) tr(JKf JK"[|)] + a t D(l - v) nKK T n, 

where the coefficients are given by 



(15) 



Eh Eh3 5 7 

C = , D = ; r- , Ct, = - , Oif — — . (16) 

1 — v 12(1 - v 2 ) ' 6 ' 10 V ; 

Here E is the Young modulus, v is the Poisson ratio, C is the stretching (membrane) 
stiffness of the plate, and D is the bending stiffness. The values of the two shear correction 
factors a s and a t from (fT7j|) have been determined in [77] / through a numerical treatment 
of several non-linear shell structures. We observe that the form (Tl5"l) of the strain energy 
density W can be obtained from the more general representation ([14"]) by choosing the 
following coefficients 

a\ = Cv= - — ^- h, a2 = 0, a 3 — C(l — u) — 2/a/i, a 4 = a s C(l — v) = 2fj,a s h, 
A + 2/i 

P, = Dv=^-^ fc = 0, h = D{l-v) = ?f t p 4 = a t D(l-v) = ^f, 
A + 2/i o oo 

(17) 

where \i and X are the elastic Lame moduli of the isotropic and homogeneous material. 
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3 Existence theorem for isotropic plates 



We employ the usual notations for the Lebesgue space L 2 (u) and the Sobolev space H 1 (uj), 
endowed with their usual norms || • \\l 2 (ui) an d || • Wh 1 ^) ■ We denote the set of proper 
orthogonal tensors by SO (3) and designate the set of (three-dimensional) translation 
vectors by M 3 and the set of second-order tensors by R 3x3 . The functional spaces of 
vectorial or tensorial functions will be denoted by L 2 (w,]R 3 ), H 1 ^, IR 3 ), and respectively 
L 2 (u>, IR 3x3 ), H 1 (u,M 3x3 ). For tensorial functions with range in £0(3), we employ the 
notations L 2 (u, SO(3)) and H 1 ^, SO(3)). We also use the classical notations for the 
norms ||v|| = (v ■ v) 1 / 2 , V« G IR 3 , and ||X|| 2 = tr(XX T ),VX G M 3x3 . 
Let us define the admissible set A by 

A={{y,Q)eH 1 {u,R 3 )xH 1 (u>,SO{3))\ y, = y*, Q, = Q*}. (18) 

The boundary conditions in ([TBI are to be understood in the sense of traces. We assume 
the existence of a function A(u, Q) representing the potential of the external surface loads 
/, c, and boundary loads n*, m* jTT] . 

Consider the two-field minimization problem associated to the deformation of elastic 
plates: find the pair (y, Q) G A which realizes the minimum of the functional 

I(y, Q) = I W{E, K) da; - A(«, Q) for (y, Q) G A. (19) 

J UJ 

Here the strain tensor E and the bending tensor K are expressed in terms of (y, Q) by 
relations (jSJ) and The variational principle of total potential energy relative to the 
functional (fl9~l) has been presented in [T7j, Sect. 2. 
The external loading potential A(u, Q) is decomposed additively 

A(u, Q) = A w (w, Q) + A duJf (u, Q), (20) 

where A w (u, Q) is the potential of the external surface loads /, c, while Ag aj/ (n, Q) is the 
potential of the external boundary loads n*, m*, which are taken in the form 



A w (u,Q) = I f -udu + U^Q), 

Jul 

Adcu f (u,Q)= / n* ■ uds + U dUf (Q). 

J du)r 



(21) 



The load potential functions n w : L 2 (u,SO(3)) -> R and U du}f : L 2 {uo, SO{3)) -> E 
are assumed to be continuous and bounded operators, whose expressions are not given 
explicitly. 

We are now able to present the main existence result concerning the deformation of 
isotropic elastic plates. 
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Theorem 2. Assume that the external loads and the boundary data satisfy the regularity 
conditions 

f G L 2 (oo, M 3 ), n* G L 2 (du f ,R 3 ), y* G H 1 {uj, M 3 ), Q* G H^oo, SO(3)). (22) 

Consider the minimization problem (|18[) . ffTU]) /or isotropic plates, i.e. when the strain 
energy density W is given by the relations (JHJ). If the constitutive coefficients satisfy the 
conditions 

2a\ + a 2 + a 3 > 0, a 2 + a 3 > 0, a 3 — a 2 > 0, a 4 > 0, , , 

2/?i + /3 2 + /3 3 > 0, /3 2 + /3 3 >0, & - & > 0, & > 0, 1 j 

i/ien the problem ffTgj) . ffTU]) admits at least one minimizing solution pair (y, Q) G 

Proof. To prove this assertion, we apply the direct methods of the calculus of variations. 
First, we observe that for any Q G 50(3) we have ||Q|| 2 = 3 and, hence, ||Q||l 2 (w) is 
bounded independent of Q. In view of the conditions f[2"2"lb 9 and the boundedness of IT^ 
and Hd idf , we derive from (1201) . (12T]) that there exist some positive constants Cj > such 
that 

|A(u,Q)| < ||/||jyi( w )||ti||ia( w ) + ||n*|| L a(a w/ )||u||La(a W/ ) + |IL(Q)| + \U dWf (Q)\ 
< Ci 1 1 l^w) + C 2 \\u\\ H i^) + C 3 + C 4 , 

which means that there exists a constant C > with 

|A(«,Q)| < C (||i/||ffi (w ) + 1), \/(y,Q)eA. (24) 

From pijl g we observe that W m b(E) is a quadratic form in the strain variables E ia (i = 
1,2, 3; a = 1,2) given by (flOj) . More precisely, we may write 

W mh (E) = W stI (E ia ) = \{a x + a 2 + « 3 )(^ 1 2 1 + E\ 2 ) + \ a 3 {E\ 2 + E\ x ) . . 

+\a i {El 1 + El 2 ) + a l E ll E 22 + a 2 E 12 E 21 . 1 J 

The quadratic form W s t r (E ia ) given by (1251) is positive definite if and only if the conditions 
j23]) 1-4 on the coefficients are satisfied. Then, by virtue of the relations ( l23]) i_/i we 
infer that there exists a constant c\ > such that 



W stT (E iQ ) > Cl ^^^, V E ia G 



or, equivalently, 

|2 



V^ m b(^) > cill^f, VE = E ia e t ®e a , E ta eR. (26) 
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Analogously, from the conditions (]2"3"]U_« on the coefficients fa we deduce that there exists 
a constant c.\ > such that 

W hend {K) > CiWKf, VK = K ia e i ®e a , K ia eR. (27) 

On the other hand, in view of (jSj) we observe that 

||.E|| 2 = tr{EE T ) = tr[(e Q ■ e p )e a <g> e p \ = s a ■ s a = y a ■ y a - 2y a -d a + 2, 

since ||d a || = 1. Then, the Cauchy-Schwarz inequality yields 

E\\ 2 du > H^illia^) + ||y, 2 |li2 H -2v / 2a(||i/ il ||i2 M + \\y, 2 \\ 2 L 2 {uj) ) 1/2 +2a 1 
or 

ll^ll£» > II^II^H " 2^ \\F\\ L2M + 2a, (28) 

where a = area(a>) and F = Grad s y = y a <g> e a is the surface gradient of deformation 
introduced previously. 

We show now that the functional I(y, Q) is bounded from below over A. Indeed, for any 
(y> Q) & A we use consecutively the inequalities fT2"4|) . (|26i) . (128]) to write 

I(y,Q)> [w mh (E)du-A(u,Q)> I Cl \\Efo\uj-C{\\y\\ m(u]) + l) 
> c i(\\ F \\h( u ) - 2 ^ ll-F|U» + 2a) - C( \\y\\HHu,) + 1) 

so that 

I(y,Q)> hWFf^-hWyWni^-h, V(y,Q)eA, (29) 

for some constants fci, fe, &3 with &i > 0, k 2 > 0. Using the Poincare inequality for the 
field y — y* 6 H 1 ^, M 3 ) (with — = on duj) in the form 

|| Grad s (y - 2/*)||| 2(a)) = || (y - y*),i||!2 H + || (y - y*), 2 ||ia( w ) > 411 V ~ v'Wmw > 

then from ( 12 9 p we deduce that there exist some constants K% > 0, K 2 > and X 3 such 
that 

I(y,Q)> ^Wy-y^l^-K^y-y^Hi^-Ks, V(y,Q)eA (30) 

Relation (]3"U]) shows that / is bounded from below over A, and thus there exists an 
infimizing sequence {y k , Q k }^ =1 C A with 

lim Q k ) = inf{ /(y, Q) | (y, Q)eA}. (31) 

fe— >oo 
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According to the hypotheses (1221) ^ we have I(y*,Q*) < oo. In view of f[3"U|) . we may 
write for I(y k , Q k ) 

oo > I(y*, Q*) > I(y k , Q k ) > K x \\y k - y*\\ 2 mH - K 2 \\y k - y*\\ HHuj) - K 3 , WkeN, 

(32) 

and hence, the sequence {y k } k=1 is bounded in H 1 ^, IR 3 ). Consequently, we can extract 
a subsequence of {y k }^_ v not relabeled, which converges weakly to an element y in 
ifV,M 3 ), i.e. 

y k ^y in H 1 ^,^), for k ->• oo, (33) 
and moreover it converges strongly in L 2 (tu,M. 3 ) by Rellich's selection principle 

y k ^y in L 2 (u,R 3 ), with y £ i/ 1 (a;,IR 3 ). (34) 

On the other hand, from ([32]) and ([19]), §Mk it follows that the sequence f u W mh (E k )du 
is bounded independent of k £ N. Then, from f l2"B"j) we deduce that {-E fc }^ 1 is a bounded 
sequence in L 2 (u, M 3x3 ). Therefore, there exists a subsequence (not relabeled) and an 
element E £ £ 2 (u;, M 3x3 ) such that 

E k ^E in £ 2 (u;,M 3x3 ), with £ fc = Q fc ' T [(y fc a - Q fc e a ) ® ej . (35) 

Similarly, from f[32]) we deduce that f Wbend(K k )du is a bounded sequence and, in view 
of (|27|) . the sequence jlif*}?^ is bounded in L 2 {u, M 3x3 ). Taking into account (EJ) and 
(EJ), we observe that ||K|| 2 = J (||Q,il| 2 + IIQ,2|| 2 ) and \\Q\\ 2 = 3- Then, it follows that 
the sequences are bounded in L 2 (u, IR 3x3 ) and, hence, {Q k }T_ 1 is a bounded 

sequence in H 1 ^, M 3x3 ). Consequently, there exists a subsequence (not relabeled) and 
an element Q £ i3" 1 (w, M 3x3 ) such that 

Q k ^Q in /f 1 (w,M 3x3 ), and Q k ^ Q in £ 2 (u;,M 3x3 ). (36) 

Moreover, we observe that L 2 (uj, SO (3)) is a closed subset of L 2 (uj, M 3x3 ). Indeed, for any 
sequence {-R fc }f\. C L 2 (u, SO (3)) which converges to an element R in L 2 (u : M 3x3 ), one 
can show that R k R T -> 11 in £ 2 (u;, M 3x3 ), and also R k R T ^ RR T in L 1 ^, M 3x3 ). It 
follows that RR T = 11 holds, which means that #£ L 2 (o;,50(3)). 
Consequently, from flUJh and Q fe £ L 2 (u, SO(3)) we obtain that Q e L 2 (u, SO(3)). 
By virtue of the boundedness of {K k }°? in L 2 (u, M 3x3 ), there exists a subsequence (not 

relabeled) and an element K £ L 2 (u, R 3x3 ) such that 

K k ^K in L 2 {u,R 3x3 ), with K k = Q k ' T [&x\{Q k a Q k > T ) ® ej. (37) 

Concerning the (weak) limits Q) -E 1 and X specified by relations fl33|) - fl37|) . it remains 
to show that they satisfy the equations 

E = Q T [(y^-Qe a )®e a \, K = Q T [ axl(Q a Q T ) <g> e a ] . (38) 
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Indeed, for any test function cp G C^°(u;,IR 3 ), we may write 

/ {Q k ' T y% - Q T y, a ) ■ = [ [{Q k ' T - Q T )v k a + Q T {y k a - y, a )] ■ 
= I [(Q k -Q,y k a ®v) + {y k a -y, a )-Q<p}duj 

Jlu 

< \\Q k - Qh*(u) \\y k a ® ¥>IU"( W ) + / (i/J, - y, a ) ■ Q^^LU. 

Taking into account fl33|) . fl36l) 9 and the boundedness of {?/ fc } in H 1 ^, M 3 ), we deduce 
that - Q||l2(w) 0, and \\y k a ® <£>||l 2 h is bounded, and y fc a — ^ y Q in L 2 (a>,IR 3 ). 
Then, from the last inequality we obtain 

f (Q k < T y k a ) ■ Vdu — ► / • ^dcu, V V G C °>,R 3 ). (39) 

But from (|35p we see that the sequence {Q k,T y k a \T_ l admits a weak limit £ in L 2 (u, M 3 ). 

By virtue of (139]) we infer that this weak limit ^ must coincide with Q y . Thus, we 
have proved the convergence 



Q k ' T y k a ^Q T y >a in L 2 (u,R 3 ), a = 1,2, 
or, equivalently, 

Q k ' T [(y k a - Q k e*) ® ej - Q T [(y, a - <3 e «) ® e «] in ^ 2 (^,m 3 ), (40) 



which means (in view of (1351) ) that the relation (j38]) i holds. 

Next, in order to prove the relation (1381) ? we proceed analogously: for any test function 
4> G C[5°(u;, M 3x3 ), we have 

= / '(Q k - Q , 0Q^) du + I {Q k f -Q T a ,Q T 4>)^ 

J U) J U) 

<\\Q k - H0Q fc a IU» + / <Q fc a - Q, Q , </> T Q > da;. 

Since the sequence {Q^ Q }^_ 1 is bounded in L 2 (u, IR 3x3 ), in view of (155]) we derive that 
the right-hand side of the last inequality tends to zero for k — > oo. Thus, we obtain 

f(Q k Q k ^,cf>)duj^ [(Q^^^du, V^C»( W ,l 3x3 ), a = 1,2. (41) 



11 



By virtue of the relations (jSj) and ([37]) ? we can write 

2 

ll^l&M = Haxl(Q fc a Q fc ' T ) (S e Q ||| 2(w) = £ || axl^Q^Hi^ 

a=l 

! 2 



2 

a=l 



Since {K k }™ is bounded in L 2 (u : IR 3x3 ), we deduce from the above relations that the 
sequences {Q k Q k ( £ r } C ^_ 1 are also bounded (a = 1,2) and, hence, they admit weak limits 

A A A X" 

in L 2 (w,IR 3><3 ). Now, in view of flUJ) we obtain that £ a = QQ a , i.e. 



Q k Q k X ^ QQ, a m L 2 (u,R^), 

or 

axl(Q fc a Q fc ' T ) ® e a axl(Q a Q T ) <g> e a in £ 2 (u;, M 3x3 ). (42) 

Taking f )36|) 9 , ( 157)) and ( I42p into account, and repeating the above argument (with arbi- 
trary test functions </> 6 C™{u), IR 3x3 ), similar to ( I4T)) ) it follows that 



Q k ' T [&xl{Q k a Q k ' T ) ®e Q ] Q [axl(Q Q ) <g> e J in L 2 (u 



,M 3x3 ), 



which means that the relation fl38l) o holds. 

In the next step of the proof, we show the convexity of the strain energy density function. 
By virtue of the conditions on the constitutive coefficients f[2"5)) . the Hessian matrix of the 
quadratic form W m h(E) in (125)) is positive definite. Similar arguments hold for Wh en d(K), 
and altogether we obtain 

W(E, K) is convex in (E,K). (43) 

By (J35]),([3ZD and ([38D we deduce E k E and K k ->> K in L 2 (w), and from (03]) we find 
that 

A") dw < lim inf / H/(£ fc , K fc ) du. (44) 

If we denote by u k = y k — x , ii = y — x, then from (l2Til) (I2"2l.(l511).(15fil).) and the continuity 
of , Hg Uf it follows 

lim A(u k ,Q k ) = A(u,Q). (45) 

k— >oo 

Since the pairs (y k , Q k ) satisfy the boundary conditions on dud , we deduce in view of the 
convergence relations (155)) . (154]) . (I5S)) and the compact embedding in the sense of traces, 
that y = y*, Q = Q* on dojd . Hence, we have (y, Q) G A. 

Finally, from (JSJ,(@U and (05]) we obtain I(y, Q) < M{I(y,Q) \ (y,Q) G .4}, which 
means that (y, Q) is a minimizer of the functional I over A. The proof is complete. □ 
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Remark 3. We observe that Q G L°°(u, SO(3)) since \\Q\\ 2 = 3, but the rotation Q 
may fail to be continuous, according to the limit case of the Sobolev embedding. Also, 
the solution deformation y G i/ 1 (o;,]R 3 ) may fail to be continuous. This is indeed an 
advantage of our formulation, since the results are applicable to a large class of shells 
(such as e.g., shells with singularity lines), and the minimizing solutions are not restricted 
by too strong regularity conditions. 

Remark 4. We notice that the conditions on the constitutive coefficients ([23]) are satisfied 
for the particular model of isotropic plates presented in Remark [JJ Indeed, taking into 
account the identification (TT7|) . we find that the inequalities (123"]) reduce to 

li > 0, 2(j, + 3A > 0. 

These conditions are satisfied in view of the positive definiteness of the three-dimensional 
quadratic elastic strain energy density for isotropic materials. Thus, the existence result 
given by Theorem [H applies to the particular plate model presented in fTfl] [77] /. 

4 Generalization of existence result 

In this section we present some variants and generalizations of Theorem [2J 
We observe that the boundary conditions imposed on the rotation Q can be relaxed or 
even omitted in the definition of the admissible set ([IS]) . For a discussion of some possible 
alternative boundary conditions for the rotation field Q on dua we refer to the works 
[27] 125] . In this line of thought, we present next the existence result corresponding to a 
larger admissible set. 

Theorem 5. Consider the minimization problem ([TO]) , over the admissible set 

A= {(y,Q) G H\u,R 3 ) x H\u,SO(3)) \ V\ dud = V*}- (46) 

If the external loads f , n* and the boundary data y* satisfy the conditions (]22]) i_^ and the 
constitutive coefficients , 0^ verify the inequalities ()23l) . then the minimization problem 
(jl9|) . (]46[) admits at least one minimizing solution pair (y, Q) G A. 

Proof. The proof can be achieved in a similar manner as the proof of Theorem [2] where 
the boundary condition Q = Q* on dud has not played an important role. □ 

The Theorem [2] is concerned with isotropic plates for which the strain energy density 
W(E, K) is given by relations ()14p . We can generalize this existence result to the case of 
anisotropic non-linear plates, provided the function W satisfies the conditions of convexity 
and coercivity: 
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Theorem 6. (Anisotropic plates) Consider the minimization problem (fTgj) . ffTPl) as- 
sociated to the deformation of anisotropic plates, and assume that the external loads and 
boundary data satisfy the conditions (|22i) . Assume that the strain energy density W(E, K) 
is an arbitrary quadratic convex function in (E,K), and moreover W is coercive, in the 
sense that 

W(E,K)>k(\\E\\ 2 + \\K\\ 2 ) 1 \f E = E ia ei®e a , K = K ia ei®e ai E ia ,K ia eR, (47) 

for some constant k > 0. Then, the minimization problem (|T8|) . (fT9l) admits at least one 
minimizing solution pair (y,Q) G A. 

Proof. We follow the same steps as in the proof of Theorem |2j In view of (T4T|) , we show 
first that the estimate ( 130|) holds also in our case. Then, there exists a minimizing sequence 
{y k , Q k }^L 1 C A, and we can prove similarly that it verifies the relations (l3T|) - (15g|) . By 
virtue of our hypothesis, W is convex in (E,K), so that the properties ( l4"5j) - (T4"o]) are 
satisfied, and we can reach the conclusion of the theorem. □ 

The Theorem [6] remains valid also for the minimization problem written over the larger 
admissible set A given by (146"j) , instead of the admissible set f fTHj) . 

Remark 7. The 6-parametric theory of shells can be used to model also composite thin 
elastic structures. In this case, the internal energy density has a more complicated struc- 
ture, and exhibits multiplicative coupling of the strain tensor E with the bending tensor 
K, see e.g. Iffi \2^. Nevertheless, the Theorem^ can be applied to deduce the existence 
of minimizers for layered composite plates, under appropriate conditions on the mate- 
rial/geometrical parameters |2J/. 

5 Comparison with a Cosserat model for plates 

In this section we present a comparison with the Cosserat model for plates (shells) pro- 
posed and investigated by the second author in [27J 129] • This model is obtained by a 
consistent formal dimensional reduction of a finite-strain three-dimensional Cosserat (mi- 
cropolar) model to the two-dimensional situation of thin plates. 

Apart from the differences in notations, there are essential similarities between this 
Cosserat plate model and the 6-parametric plate theory presented in Section [TJ Firstly, 
in both approaches the primary independent kinematical variables are the deformation 
field y G M 3 and the rotation tensor field Q G SO(3), which are denoted in by m 

and -R. Thus, we have the correspondence 

y = m, Q = R. 
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The rotation field R in this derivation approach is inherited from the parent Cosserat bulk 
model, which already includes a triad of rigidly rotating directors (the "triedre mobile", 
see [33]). 

Moreover, the measures of strain are essentially the same in the two approaches. Indeed, 
the so-called stretch tensor U and the third order curvature tensor & s are introduced in 
[2"T] Sect. 4, through the relations 

U = R T (m^ a (g> e a + R 3 <g> e 3 ), 

& s = [R (Ri, a <8> e a )l <g> e-i , with i?j = Rei . 

If we compare these tensors with the strain measures E and K defined in (J8])- (|lip . then 
we find that 

E = U-TL, 

where II = <g) is the identity tensor. It can be shown that the tensor ^ s has only 6 
distinct components which coincide with the 6 non-zero components of the bending tensor 
K given by ( ITT]) . 

Let us compare now the expressions for the strain energy density. The strain energy 
density is assumed in the Cosserat plate model in the form 

W(U,ft s ) =h(ji ||sym(C7 - 11) || 2 + ^||skew(t7 - 11) || 2 + — ^ (tr(C7 - E)) 

+ ^ (/i||sym^||2 + Mc || ske w^|| 2 + (tr(i^)) 2 ) (4J 



+ [ 2 (l + «4^||^ s || 9 )^(«5||sym^|| 2 + a 6 ||skew^|| 2 + a 7 (tr(^)) 

i=i 

where 8l s = & s ei , and fj, c > is the Cosserat couple modulus, L c is an internal length, 
and a 4 , ...,a 7 are constitutive coefficients introduced in [27]. The exponents p and q are 
such that p > 1, q > 0. The existence of minimizers is proved for the case of the strain 
energy (j4"8|) under the simple conditions 

/x c > 0, /i > 0, A > 0, a 5 > 0, a 6 > 0, a 7 > 0. (49) 

If we choose the exponent p = 1 and the parameter = in (|48i) . then the expressions 
of the strain energy density in the two approaches coincide, in terms of the independent 
kinematical variables (y, Q) = (m, R). To realize the coincidence of the two strain energy 
functions (j!4p and 0481) in this case, we need to identify the set of constitutive coefficients 
(«i, a 4 , /3 4 ) from 0T41) with the set of parameters (A, fi, /i c , L c , a 5 , a 6 , a 7 ) in the 



3 i+p 
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following way 

ax = h 2 ^ , a 2 = h(fi - fie) , a 3 = h(fi + /i c ) , a>4, = K>h(fi + /j, c ) , 

A + Zfi 

ft = P 1 *^ - ^) + ^ - , ^ = ~ ( h2 y^ + L l a , (50) 
ft = -y [^^^7 + (| «5 + | a 6 + a 7 )] , ft = ^ L\ (§ a 5 + | a 6 + a 7 ), 
where fc is the formal shear correction factor. We observe that 

a 3 - a 2 = 2hfi c , (51) 

so that the condition \i c > considered in (149"]^ corresponds to a 3 — a 2 > 0, assumed in the 
conditions f[23]) which ensure the positive definiteness of W. The interesting degenerate 
case fi c = is investigated in details in [29]. In view of fl5T|) . this case corresponds to 
«3 — ^2 — 0, when the energy function W is only positive semi-definite, and the proof of 
the existence results is more delicate [29J. 

There are also some differences between the two approaches, such as for instance the 
form of the boundary conditions and the expression of the external loading potential 
A(u,Q). The proof of Theorem |2] follows the same steps as the proof of Theorem 4.1 
in [27]. However, the conditions on the constitutive coefficients ( 149]) imposed in [27] are 
more restrictive than the conditions ( 12 3 p in Theorem HJ For example, in the case of the 
isotropic plate model described in Remark [[] the constitutive coefficients ([17]) satisfy the 
existence conditions ( 123]) stated by Theorem [2] but they do not verify the restrictions (T49]) 
(in virtue of the identification (150]) ). This illustrates the fact that the conditions for the 
constitutive coefficients established in the present work are less restrictive. 
Finally, we remark that the form of the strain energy density W given by (145]) is more 
general than f[T4]) . Indeed, the expression ( 1141) of W can be obtained from ( I48p if we 
choose the parameters p = 1 and 04 = 0. Nevertheless, the existence of minimizers is 
proved in [27] for any exponents p, q and coefficient 04 satisfying the conditions: p > 1, 
q > 0, and 04 > 0. 

In a future contribution we will extend our results to the general case of 6-parameter 
shells. 
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